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Abstract 



We study one-dimensional reaction-diffusion models described by 
master equations and their associated two-state quantum Hamiltoni- 
ans. By choosing appropriate rates, the equations of motion decouple 
into certain subsets. We solve the first subset which has a close rela- 
tion to the problem of lattice electrons in an electric field. In this way 
we obtain L(L — 1) + 1 energy levels of a quantum chain with L sites. 
The corresponding Hamiltonian depends on 7 parameters and does 
not look integrable using conventional methods. As an application, 
we compute the dynamical critical exponent of a new type of kinetic 
Ising model. 



1 Introduction 



One-dimensional quantum chain Hamiltonian show up in different con- 
texts in equilibrium statistical mechanics and many examples of integrable 
ones are known. They appear as the logarithmic derivative of the mon- 
odromy matrix related to the Yang-Baxter equation [0J. However, there is 
also a connection to stochastic problems, because a master equation can be 
related to a suitable Schrodinger equation ||. For example, the diffusion of 
classical particles with hard-core repulsion on a lattice leads to the ferromag- 
netic Heisenberg model ||, 0. It was realised recently ||, [|], 0, |§ that in 
various other stochastic models the corresponding Hamiltonian is integrable 
and one can use this integrability in order to determine the phase diagram of 
the system and compute the time evolution of various average quantities. In 
the Schrodinger equation the wave function is the time dependent probabil- 
ity distribution, the Hamiltonian is in general non-hermitian and its matrix 
elements are constrained by the condition that the sum of rates for a given 
process have to add to one (conservation of probabilities). 

In the present paper we adopt the opposite strategy. We ask ourselves if 
there are not known examples in non-equilibrium statistical mechanics where 
one can solve the master equation or part of it (in a sense to be defined soon) 
and thus get in this way the spectrum and wave functions of a Hamiltonian. 
We first looked at the Hamiltonian related to the master equation studied by 



Kimball [Bl and Deker and Haake 10 1 where two wave functions were known 



for any number of sites, a few more can be found |Tl[ but we were not able to 
go beyond this point. Another approach which was proven very powerful is to 
consider various moments of the probability destribution (the knowledge of 
all the moments giving back the probability distribution) and to find special 
rates in the master equation such that the differential equations for various 
moments decouple into subsets. The first example of this kind is the Glauber 
model [12 1 used to study the relaxation of the one-dimensional classical Ising 
model. 

Going to the dual variables description ||13|| , the Glauber model corre- 
sponds to solve the one-spin-string problem (see Sec. 2) which decouples 
from the many-spin-strings problem. Subsequently [[14]], |T5j the many strings 
problem was also solved and implicitely the whole spectrum of the Hamilto- 
nian derived. The Hamiltonian is however trivial: it corresponds to the XY 
model in a Z field which can be brought through a Jordan- Wigner transfor- 
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mation to a bilinear problem in fermionic creation and annihilation operators 
and can thus be trivially diagonalized. 

Another case is the one considered in Refs. JIE| , JT7J, [0 and [[TJJ. This 
is an example where the master equation describes chemical processes and 
brownian motions [ 20(| . In this case (see Sec. 2 for the terminology) one 



considers first processes where two identical molecules A give a molecule 
A and a vacancy (coagulation) with a rate "tuned" to the diffusion rate 
and the reverse process (decogulation) with an arbitrary rate (A+ vacancy 
— > A + A). This is a one-parameter process since one rate can be used to 
fix the time scale. In this case one considers a one-hole probability function, 
two- holes etc. ... as moments of the probability function. It was shown that 
the equations for the one-hole probability function decouple from the others. 
The Hamiltonian in this case is again that of the XY model in a Z field. 
Some work was done also for the two- holes problem [BTJ , ||22||. Actually, as 



will be shown in this paper knowing the two spin-string solution of |T3| and 



y| for the Glauber problem solves the two-holes problem of the coagulation- 



decoagulation case. The Glauber model corresponds to chemical processes 
in which two molecules annihilate into two vacancies: A + A — > vacancy + 
vacancy with a rate tuned to the diffusion rate and the reverse process with 



an arbitrary rate. The investigation of Refs. |16| - |19| has, however, one 
more case, namely to the coagulation-decoagulation processes one adds the 
birth process (two vacancies give a molecule A and a vacancy). They were 
able to solve the continuum limit of this problem in the case of one-hole. The 
spectrum can be expressed in terms of zeroes of one Airy function and the 
Hamiltonian is not any more a trivial one. This observation triggered the 
present paper which is organized as follows. 

In Sec. 2 we give the master equation for a chain of L sites, having on 
each site a stochastic variable taking two values. The rates are assumed 
to depend on the configuration of two neighbouring sites only. Next we 
present the dictionary necessary for chemistry which is also useful in order 
to understand our results. The one-dimensional quantum chain associated 
with the master equation is also given. 

In Sec. 3 we first define the holes and spin-strings. Next we derive 
the conditions under which the differential equations for one-hole decouple 
from the others. These conditions assure also that the two-holes problem is 
decoupled from the three or more holes, etc. since there are five conditions for 
the twelve rates. This leaves the Hamiltonian dependent on six parameters 
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(the normalisation fixes another one). The Hamiltonian one obtains is non- 
hermitian. We write the differential equations for the closed ring case and 
notice that they depend on one parameter less. This is not the case for 
the open chain. We also show that the same differential equations describe 
the one-string problem and give the relation between the parameters which 
appear in the differential equations and the rates of the one-string problem 
(the conditions on the rates are different in this case but their number is of 
course the same). 

In Sec. 4 we solve the one-hole subset of differential equations for a finite 
chain on a closed ring. Their structure depends on one parameter (called 5 in 
the text) which is non-negative in the stochastic problems. If 5 is non-zero, 
the eigenvalues are given by zeroes of Lommel polynomials. If 5 is zero, the 
eigenvalues are of trigonometric type. We also check under which conditions 
the spectrum corresponds to a free fermionic one. 

The continuum limit is considered in Sec. 5. If 5 ^ one can take the 
limit L —>■ oo and find that like in the special case of Refs. [l6j - Jl9| the 



eigenvalues are related to zeroes of one Airy function. For 5 = 0, they have 
a simple expression. The eigenvalues are in general complex. 

As a simple application of our approach, we consider in Sec. 6 the critical 
dynamics of the Ising model [[y|, pjjl , [f24| . This process is described by 



the same master equation as the one used for diffusion-reaction processes 
with supplementary conditions on the rates which assure that for large times 
the system reaches the equilibrium distribution of the one- dimensional Ising 
model (detailed balance conditions). We have checked whether our seven 
parameter Hamiltonian is compatible with the detailed balance conditions 
and found two solutions: one is the old Glauber one | i~2|| , | I3]| , another one 
is new. Its physical signifiance is discussed. 

One last comment before going to the bulk of this paper: one-dimensional 
diffussion-reaction processes are apparently seen experimentally This 
observation might make the reader have a closer look at the subject. 
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2 The Master Equation and 1 - d Quantum 
Chains 



We consider a one-dimensional chain with L sites and on each site k(k = 
1, 2, ...L) we take a variable f3 k which takes two values: and 1. For f3 k = 
the site k is empty (vacancy), for (5 k = 1 the site k is occupied by a molecule 
A. At t — the probability to find a certain configuration of molecules and 
vacancies on the chain is given by the probability distribution 

P o (fr,&,...ft)=Po({0}) (2-1) 

If we assume that we have two-body processes only, the time evolution of the 
system is given by a master equation: 

J2l-V((3 k ,(3 k+1 )P((3 u ...,(3 L ;t) 
k=i 

i 

Y! r£3£i p (&» • • • » A-i, 7*. 7fc+i, • • • , ft; 01 

7fc,7fe+l=0 

(2.2) 

which gives the probability distribution at the time t for a given P({fi},t = 
0) = Pq. The first term on the r. h. s. of (2.2) describes the losses and the 
second one the gains. A prime is used to indicate that in the sum the couple 
Ik = Ac, 7fc+i = Pk+i is excluded. T^'^^ represents the probability per unit 
time that the configuration (7fc,7fc+i) on neighbouring sites changes into the 
configuration (/3 k ,f3 k+ i). The conservation of probabilities gives: 

v(i k ,i k+ i) = £' r-;?- (2-3) 

A>/3fc+i=0 

where again the prime in the sum indicates that the couple ^ k = f3 kl ^ k+ i = 
Pk+i is excluded. The master equation (2.2) is quite general and can have 
various physical interpretations (in Sec. 6 we will apply is to study critical 
dynamics), in the present section however we use the language of diffusion- 
reactions processes where the rates have a simple meaning. The following 
processes are included in the master equation (with a vacancy (f3 = 0) and 



dP (/¥) 
dt 



+ 
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A a molecule ((3 = 1)) 



(2.4) 



Diffusion to the right: A + -> + A (rate rj?) 

Diffusion to the left: + A -> A + (r?J) 

Coagulation at the right: A + A — > + A (T^\) 

Coagulation at the left: A + A -> A + (T\l) 

Decoagulation at the right: A + A + A (T\^) 

Decoagulation at the left: + A -> A + A (T°\) 

Birth at the right: + -> + A (r°°) 

Birth at the left: + -> A + (r?°) 

Death at the right: + A -> + (r{£) 

Death at the left: A + -> + (rjg) 

Pair-annihilation: A + A -> + (rJJ) 

Pair-creation: + -> A + A (rj?) 

In Ref. (see also the references included therein), several special cases 
have been considered in detail. The notations for the rates used in Ref. 
are related to those used in the present paper: 

rjs = Wj-aj-pia, P) ; V(a,P) = w 00 {a,(3) (2.5) 
To a given master equation one can associate [0] a Schrodinger equation 

d_ 

Of 

as follows. Take a 2 L -dimensional orthogonal basis: 

| {(3} >=\ p 1 ,p 2 ...,f3 L > , < {/?'} | {(3} >= 8 w>m (2.7) 

and denote: 

| P > = £ P({/3}; t) | {/?} > ; | P > = E P o({/9}) | {/?} > (2.8) 

P P 

In order to write the Hamiltonian, we first define a basis in the space of 2 x 2 
matrices taking E a '^(a,f3 = 0, 1) with matrix elements: 

(E a ^ s = 5 ar/ S PtS (2.9) 



P > = -H\P> (2.6) 
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On each site k of the chain we take the matrices E^ ,l3 (k = 1,2, ... L). It was 
shown in Ref. that in this basis the Hamiltonian H occuring in Eq. (2.6) 
has the following expression: 



H 



L-l 
k=l 



where 
with 



pa/3 
" 7,0 

r a,/3 
a,f3 



E7**,/8 7717" J?&P 

7, (5 = 5 excluded) 



-V(a,(3) 



(2.10) 

(2.11) 

(2.12) 
(2.13) 



where we have used Eq. (2.3). The Hamiltonian H is in general non- 
hermitian and due to the probability conservation (Eq. (2.3)), it satisfies 
the remarkable relation: 

<0 \H = (2.14) 
where the bra ground-state < | is 



< 



£<{/?} 





(2.15) 



From Eq. (2.14) it follows that the ground-state energy is zero. The Hamil- 
tonian might have more than one eigenstate of zero energy i. e. there may 
be more than one stationary solution. 

It is instructive to write the 4x4 matrix given by Eq. (2.11) explicitely: 



/ roo i poo i poo 

/ 01 10 11 



00 



^00 



pOO 
_i 01 

pOO 
1 10 

pOO 

_i 11 



pOl 
1 00 



_p01 
1 00 

pOl 
1 10 

pOl 

1 11 



pOl 

1 11 



plO 
1 00 



_pl0 
1 00 

plO 
1 01 

-r ?H 

plO 

1 11 



plO 

1 11 



_pll 

1 00 

_pll 

1 01 

_pll 

1 10 

pll I pll I 

1 00 ^ 1 01 ^ 



pll 

1 10 



V 



(2-16) 

Notice that the sum of each column in (2.16) vanishes, this is a consequence 
ofEq. (2.14). 



/ 
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Let us now consider an observable X (fti, . . ./?£,) = X({(3}), its average 
quantity can be computed as follows: 

<x>(t) = £x({/3})P({/?},t) 

= < I X I P >=< I Xe- Ht | P > (2.17) 
The time derivative of < X > (t) is obviously: 

- < X > (t) = - < | XHe~ m | P > = < | [iy,X]e- Hi | P > (2.18) 

where we have used Eq. (2.14). In the next Section two types of observables 
will be used. The first one is a hole of length n starting at the site j: 

j+n-l 

X= II b{{3 k ) (2.19) 
k=j 

Its average is: 

K(j,n + j-l) = Y, S ([ 3 jWj+i)---t(Pj+n-i)P{{f3};t) 

n+j-l 

= < I [] E° k °e- m | P > (2.20) 



k=j 



site j, 



where we have used Eq. (2.9). The time dependent function K(j,n + 
1) gives the probability that n consecutive sites starting with the site j, 
are empty (a hole). This observable is sometimes used also in equilibrium 
problems p6| . A second observable which will be of interest is a spin-string 
of length n starting at the site j: 



j+n-l 



X= II W k )-6(l3 k + l)) (2.21) 

k=j 



Its average quantity is: 



n+j-l 

S(j, n + j-l) =<0| n °l e ' m I p o > (2-22) 

k=j 



7 



where 



= E 00 - E 11 (2.23) 



Actually one can define an observable interpolating between the one-hole and 
the one-spin-string taking 

x = 3 jf {ih my ih m + 1} ) (2,24) 

where £ is a parameter. Its average quantitiy being: 

n+j-l 

X(j, n + j-l) =<0 | II ^ I P o > (2.25) 

k=j 

where 

r = ^E 00 - -^E^ (2.26) 

3 Decoupling of the master equation into sub- 
sets 

As we have seen in the proceeding Section, the master equation (2.2) is 
equivalent to the Schrodinger equation (2.6). Thus finding the solution of 
the 2 L linear differential equations (2.2) will give the eigenvalues and eigen- 
functions of the Hamiltonian given by the Eqs. (2.10) - (2.11). The main 
idea that we are going to use in this paper is to find conditions under which 
the 2 L differential equations decouple into subsets which can be solved inde- 
pendently. This corresponds to bringing the Hamiltonian to a block-diagonal 
form. There are some obvious sub-sets that one might consider: 

a) Holes The definition of a hole was given by Eq. (2.20). One can ask 
for conditions on the rates such that the one-hole differential equa- 
tions decouple from the others. The two-holes equations decouple from 
more than two- holes etc.... Since obviously knowing the solution of all 
the holes probabilities gives the probability distribution P({(3},t), this 
corresponds to the diagonalisation of the Hamiltonian. 
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b) Spin-strings The definition of a spin-string was given by Eq. (2.22). 
One can consider the hierarchy of the sub-sets of one string, two strings 
etc. Again knowing all the spin-strings functions is equivalent to the 
knowledge of P ({/?}, t) (the proof is trivial). 

c) Vacancies The iV-vacancies function is 

V(h,i 2 ,...i N ) = £5(AJ5(AJ...5(A JV )P({/3};t) 

P 

= <0|«...< e- Ht \P > (3.1) 



The hierarchy is one vacancy, two vacancies etc. ... (see also |27[| , [< 
11)- 



d) Spins The iV-spins function is 

p 

= <0|<<...< JV -«|P > (3.2) 

with the hierarchy one spin, two spins etc. ... 

As we are going to show soon, once we solve the holes case we can get through 
a similarity transformation the spin-strings solution. 

It is interesting at this point to solve a small counting problem. Take first 
an open chain with L sites and count how many kind of holes N m (various 
lengths and positions) one can have for a given total of m holes (m — 1,2 
for L even, m — 1, 2, . . . ^±1 for L odd). The answer is: 

N m = + = C 2 L m + C 2 /"- 1 (3.3) 

(2m)!(L-2m + l)! L L K ' 

The meaning of the right hand side of Eq. (3.3) will be apparent imme- 
diately. Consider now for example the following Hamiltonian: 

H = L E H t (3.4.a) 

i=l 



L 
' •' 2 



H = - 2 



Vvfvf+i + V X °\°\ + °t + 1 -V-l] (3-4.b) 
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where 77 is a free parameter and a x , a y and a z are Pauli matrices. This Hamil- 
tonian will correspond to a special choice of the parameters in our solution 
of the decoupling problem. Here we are only interested in the structure of 
the problem. The Hamiltonian (3.4) can be diagonalized through a Jordan- 
Wigner transformation [3Cfl. The result for the open chain has the form: 



L 

H = J2 £ kn k (3.5) 

k=l 

where n& are fermionic occupation numbers (n& = 0, 1) and E\ = (there is 
a zero fermionic mode). If one computes now the number of states in the 
sectors with 2m and (2m— 1) fermions one finds precisely (3.3). The fermionic 
vacuum does not appear in the holes calculation. The reason is simple: 
in a probabilistic calculation one does not need to know 2 L probabilities 
P(Pi, • • • , Pl) but only 2 L — 1 since the sum of probabilities adds to one. 

We now consider a closed ring with L sites. In this case the num- 
ber of holes in the subset with m-holes (m — 1, 2, . . . , £ for L even, m 
1,2,..., for L odd) is 



' 2 



2L\ 

Nm = (2m)\(L-2m)\ (3 ' 6) 

except for m = 1 where we have one more hole than one obtains from Eq. 
(3.6). This special hole corresponds to the situation where all the L sites are 
empty. The reason for this choice of the supplementary hole comes from a 
decoupling of this hole from the other functions in the 1-hole subset (see Sec. 
4). Adding all states coming from all the m-holes subsets one gets 2 L — 1 as 
one should (one more state comes from probabilities conservation). 

We now return to the Hamiltonian (3.4) with periodic boundary con- 
ditions which can be diagonalized through a Jordan- Wigner transforma- 



tion |pl| , p2[ and its spectrum can be obtained taking the even number of 
fermions sector of two free-fermionic Hamiltonians, one with periodic bound- 
ary conditions and the other with antiperiodic boundary conditions. Notice 
that both Hamiltonians have a zero ground state energy. The number of 
states in the 2m fermions sector of both Hamiltonians give together N m of 
Eq. (3.6). We are left with two fermionic vacua. One corresponds to the 
total empty one-hole the other to the redundand state (the probabilities 
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add to one). This calculation suggests that for special choices of the pa- 
rameters of our problem we will get Hamiltonians which correspond to free 
fermions. This will be proven to be indeed the case. We will also find an 
unexpected phenomenon, where the closed chain problem will be described 
by free fermions but not the open chain. We will also get that in general 
the spectrum is not given by free fermions at all, but we believe that our 
counting could help clarify the algebraic structure of the problem. 

The vacancies problem will not be considered here but can be found in 
29fl . One finds a non-hermitian Hamiltonian depending on ten parameters 
whose spectrum is given by the XXZ Heisenberg chain in a Z field. The 
same similarity transformation used for the holes-spin-strings correspondence 
can be used to solve the spins problem once the vacancies solution is known. 

We now start with the holes problem and ask for the conditions on the 
rates such that the time derivative of a one hole probability function 

j+n—l 

< o I //• II 

k=j 



dK{j,j+n-l) 
dt 



-Ht 



Pn > 



can be expressed in terms of one-hole probability functions only, 
forward calculation leads to the following five conditions: 

0; 



(3.7) 
A straight- 



en 
1 oo 



-nOl 
1 00 



-nlO 
1 00 



1 10 



pOl . 
1 10 i 



plO 
1 01 



pll 

1 01 



(3.8) 

The physical meaning of the conditions (3.8) is the following. All processes 
creating A molecules (creation, birth and decoagulation) are allowed. From 
the processes which "kill" molecules only coagulation processes are allowed 
and they are tuned to the diffusion rates: the rate of coagulation at the 
right (left) is equal to the diffusion rate to the right (left). To shorten the 
notations, we will denote: 



r io 
1 01 

-.10 



r ii. 

1 01i 



r 01 
1 10 



pll 
1 10 



a R 

d R = rg, d L = r°{ , b R = r°° , b L = r°° , c = r°° (3.9) 

Thus the a's describe diffusion and coagulation, the ef s decoagulation, the 
fe's birth and c pair-creation. If we take into account Eqs. (3.9), the 4x4 
matrix Hk of Eq. (2.16) becomes: 

{ c + b L + b R 

-b R a L + d L 

-b L -a L 

-c L ~d L 



\ 





-or 
a R + d R 

-d R 





CLL + a R J 



(3.10) 
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and it is the Hamiltonian (2.10) with (3.10) that we want to diagonalize. 
This Hamiltonian is non-hermitian and depends on seven parameters, one of 
which will fix the time scale. It is now useful to use the chemical picture 
given by Eq. (2.4) in order to look for global symmetries of the Hamiltonian. 
In general one finds none. All one can say is that if the rates are all positive, 
all the eigenvalues are non-negative. 

The differential equations derived from Eq. (3.7) for the one-hole proba- 
bility functions for a chain with L sites depend on the boundary conditions: 
Closed ring, L(L — 1) + 1 equations 



dK(j,j+n-l) 
dt 



dK(j,j + L-l) 
dt 

with the notation: 



= -[ 1+ ( n -l)5}K(j,j+n-l) 

+a L K(j + l,j + n-l) 
+ a R K(J, j + n-2) + a L K(j, j + n) 

+ a R K(j-l,j + n-l) (1 <n < L - 1) (3.11) 

= -L5K(j,j + L-l) (n = L) (3.12) 



K(j,j) = l (n = 0) (3.13) 

Obviously K(j,j + L — 1) is j independent. 
The open chain, L(L + l)/2 equations. 

Equation (3.11) stays unchanged except when j = 1 or j = L — n + 1 
when we have: 



dK(l,n) 
di 



dK(j,L) 
dt 



= -[iL + (n-l)6]K(l,n) + a R K(l,n-l) 

+a L K(l,n + l) (3.14) 

= -h R + (n-l)5]K(j,L) + a L K(j + l,L) 

+a R K(j-l,L) (3.15) 



Equation (3.12) is replaced by: 

dK(l,L) 



di -(L-1)5K(1,L) (3.16) 
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and the notation (3.13) used in Eqs. (3.11), (3.14) and (3.15) stays un- 
changed. In the above Equations we have denoted: 





= a>R + d R -b R - c = 


-nlO , -nlO 
01 11 


a L 


= a L + d L -b L -c = 


-pOl i -pOl 
10 11 


a R 


— r 10 n — r 01 

— 1 01' U L — 1 10 




1R 


= a L + a R + d R = T°Q ■ 


, r io , r io 

01 11 


1L 


= a L + a R + d L = r?o ■ 


, r io , r oi 

^ 1 01 ^ 1 11 


7 


= lR + lL 
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= b R + b L + c = T™ + T™ + T™ 



pOO _ pOO 
1 11 1 01 
pOO _ pOO 
1 11 1 10 



(3.17) 

Notice that the equations for the open chain depend on seven parameters 
(a R , OiLa Rl a^, 7_r, and 5). In the case of the ring however, the differential 
equations depend only on six parameters since j R and appear in the 
combination 7_r + 7l only. One also sees that the time derivative of the pro- 
bability to find a hole of length n is related to the probability functions for 
one-hole of length n — 1, n and n + 1 . 

We now consider the two-holes probability function: 



n+j—l jr'+m— 1 

K(j, n+j-1 ; k,k+m-l) =< | f[ E° h ° ]J E° s ° e~ m \ P > (3.18) 

r=j s=k 

which describes the probability of having a hole of length n starting on the 
site j and a second hole of length m starting on the site k. A straighfor- 
ward calculation gives the following differential equations for the two-holes 
probability function in the case of the ring geometry: 

dK(j,j + n — 1; k, k + m — 1) 
di 

- [27 + (m + n - 2)5] K(j,j + n - 1; k, k + m - 1) 

+ a L [K(j + 1, j + n - 1; k, k + m - 1) + K(j,j + n - 1; k + 1, k + m - 1)] 

+ a R [K(jJ + n - 2; k, k + m - 1) + K(j,j + n - 1; k, k + m - 2)] 

+ a L [K(j,j + n; k, k + m - 1) + K(j,j + n — 1; k, k + m)] 

+ a R [K(j - l,j + n - l;k,k + m - 1) + if(j, j + n - 1; k - 1, A; + m - 1)] 



(3.19) 



13 



As one notices the conditions (3.8) which decouple the one-hole problem from 
more holes, decouples the two-holes problem from three and more holes. The 
two-holes problem contains one-hole probability functions. This happens 
when functions of the type 

j + n - 1; j + n, k + m - 1) = K(j, j + m + n - 1) (3.20) 

occur in the right hand side of Eqs. (3.19). The same pattern occurs in the 
problem of three or more holes and thus the conditions (3.8) have produced 
a decoupling into sub-sets of the master-equation. The same happens for the 
open chain. 

We now turn to the decoupling problem in the case of spin-strings. The 
one-spin string is defined by Eq. (2.22), the probability for two-spin-strings 
is: 

j+il— 1 fc+ra— 1 

S(j,j + n-l;k,k + m-l) =< | JJ a z r JJ a z s e~ Ht \ P > (3.21) 

r=j s=k 

etc. ... Let us show that the decoupling problem for the spin-strings takes 
us to the same differential equations as the ones for the holes with a simple 
replacement K — » S. The expressions of the constants a^z,, b^L, d^L and c 
respectively aR,L,lR,L and 5 in terms of rates have to be changed. 
Consider two matrices b and b^ 1 : 

- (v. 1 ) ^- Co 3 

and correspondingly the matrices: 

B = J] b k , 1 = n ^ (3-23) 

k=l k=l 



we perform a similarity transformation [[3^] which takes us from the Hamil- 
tonian H given by Eq. (2.10) with given by Eq. (3.10) to H s respectively 

H s = B~ X HB (3.24) 

and 
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Hi 



I b R +b L , 3c 
2 " r 4 s 
_ f M J- 

e 

4 



(b R +b L ) _ 3c i 
2 4 2 

+ a L + f 



2 
2 



rfz, I c 

2 4 



(bR+b L ) , 

2 



2 

bR-d R 



3c 
4 



2 



2 



2 -O-R + Z 



_dg_ , c 
2^4 



fefi I d,R — d,L c 

2 2 4 

_ bh_ _i_ ^ -rffl _ c 

2 2 4 

+ a L + a R - | 



(3.25) 

Notice that as in Eq. (2.16) the sum of the matrix elements in each column 
are zero and consequently the matrix elements of (3.25) can be interpreted as 
rates provided that they are positive, we will denote them by f^g. Instead of 
the five conditions (3.8) on the rates Tj^ we now find five different conditions 



on the rates T 



7,S . 



-nil , pOO _ 

00 11 ~~ 

plO i plO _ 

1 nn T 1 ii — 



f 01 _l r 10 

1 n ' 01 



00 

fill i flO _ f 01 i fSOO 
01 11 — 1 



10 
■pOO 



+ r 



10 

+ r 

11 



00 
01 



01 



pOl i pOl _ 
1 nn T 1 1 1 — 



00 

pll | pOl _ 
1 10 ~r 1 11 ~~ 



■pOO , pOO 
linTl 



10 
^00 
10 



01 
10 



(3.26) 



The spectrum of the Hamiltonian (3.25) is again non-negative if the T^l are 
positive. This implies a different domain for the parameters or^, & r ; l, d^L 
and c than the one obtained from the positivity of the rates occuring in 
Eq. (3.10). 

Let us now consider the one-hole function (2.19) and make use of the 
similarity transformation (3.24): 



j+n— 1 



K{j,n + j - 1) = < | [] E° k °Be- H 'B- 1 \ P > (3.27) 

k=j 



Let us now notice (see Eq. (2.15)) that 

< I b = < 



(3.28) 
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< I EfB = < | (E° k ° - El 1 ) = < | (E° k ° - El 1 ) =< | a% (3.29) 
With (3.28) and (3.29) we get instead of Eq. (3.27) 

j+n-j 

K(j, n + j-1) =<0 | n ^ e "^' I P o > (3-30) 

fc=j 

where 

| pS > = 5-1 l p Q > (3.31) 

Comparing the Eqs. (2.22) and (3.30) we learn that the one-hole probability 
function at time t given by a Hamiltonian H with an initial condition P ({/9}) 
coincides with the one-spin-string function given by the Hamiltonian H s and 
initial condition Pq ({/?}). The rates r^t and f^t are both dependent on the 
seven parameters or,l, &r,l, rfR,L and c but as we mentioned the domain of 
the parameters which makes the rates r^t positive is in general different 
from the domain where f ^ are positive. The differential equations satisfied 
by the spin-string functions which depend on the parameters or^, o^l, 1r,l 
and 5 are identical with those satisfied by the holes probability functions, 
the expression of these parameters in terms of the rates f ^ however, differs 
from (3.17) where the rates were used. For example: 

5 = 2(fSg + fjg) ; 7 = 2(fS + f 1° + f?° + fj») (3.32) 

Notice that 7 and <5 are non- negative as for the rates r^t (see Eq. (3.17)). 
Obviously the transformation (3.23) takes the two-holes probability functions 
into the two-spin-string functions etc.... The same similarity transformation 
takes vacancies (Eq. (3.1)) into spins (Eq. (3.2)). Finally let us note that 
the more general transformation 

b =('f T f) (3-33) 

takes the one-hole function into the more general observable (2.25). This 
observation is relevant for two reasons. On one side it shows that more 
general observables with the same structure bring nothing new from the 
point of view of integrability and on the other side, for each value of £ one 
obtains new rates which, when positive, give a spectrum with a non-negative 
real part. 
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4 The spectrum in the one-hole subset. Closed 
ring 

Our aim is to find the spectrum of the Hamiltonian H given by Eq. (3.10) 
in the one-hole sector. This implies solving the equations (3.11) - (3.13) in 
the case of the closed ring. For the special case of the left-right symmetric 
coagulation-decoagulation model: 

a L = a R = a; d L = d R = d; b R>L = c = (4.1) 



or 

otR = «l = a + d ; 7# = jl = 2a + d , 5 = (4.2) 

the one-hole probability equations were solved for both the ring and the open 
chain case [p^j . The spectrum is that of free fermions. One can easily show 
that in this particular case one can use a similarity transformation in order to 
bring the Hamiltonian given by (3.10) into the form (3.4) with i] = (l + d) 1//2 
(we choose a = 1). The Glauber model [12| which again is described by free 
fermions [[13] corresponds to the choice: 

dL = clr = cl , ol = or = —d ; dR = di = d , c = 2d (4.3) 



or 



ocr = a L = a , 7 R = 7 L = d , 5 = (4.4) 
As mentioned already, the Glauber model was solved in the framework of 



the one-spin-string problem. One other special case was solved in Refs. |$3 
JM, which again is a free fermionic case: 



or 



a L ^a R , b RiL = d RiL = c = (4.5) 



a L = a L , a R = a R , 7^ = 7 L = a L + a R , 5 = (4.6) 



We now consider the general case. It is convenient to make a change of 
variables: 

Ti — 1 

K(j, j + n - 1) = R(n; 3 + = R(n; p) (4.7) 
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which implies to consider the length of the hole and its center as variables (p 
is integer or half-integer). We next take the Fourier transform of R(n;p) : 



L-l 



q=0 



and get 



1 (dQ(n;q) 
D(q) \ Jt 



+ [(7 - 5) + n5]Q(n; q) = u{q)Q{n - 1; q) 



(4.8) 



+u- 1 (q)Q(n + l;q) (q = 0, 1, • • • L - 1) 

with the boundary conditions 

Q(0;q) = Q(L;q) = (q + 0) 
Q(0;0) = 1 ; Q(L; 0) = Ae~ L5t 



(4.9) 



where 



u(q) = [ 



an + Ql + i{oR - a^tg^ , 1/2 
a R + a L + i(a R - a L )tg^f 



(4.10) 
(4.11) 



D(q) = [(a R + a L )(a R + a L )] 1/2 



cos 



irq 



(a>R 


- a L \ 






\ ■ 2 










)sm — 


\a R 


+ a L J 


\CtR + 







1 . 2nq / a R - a L a R -a L 
+ - sin — — 1 

2 L \a R + a L a R + a L 



n l/2 

(4.12) 



and A is an arbitrary constant. The boundary conditions (4.10) for the 
q = case give an inhomogenous system of differential equations. Since in 
the present paper we are interested in the spectrum only, we consider only 
the homogenous equations and take 



Q hom (n; q)=u n (q)F(n,q)e- xt 
where A are the eigenvalues. The F(n; q) satisfy the equation: 
(E - n)F(n, q) = V(F(n + l;q) + F(n - 1; q)) 

where 

F(0;q)=F(L;q) = 



(4.13) 

(4.14) 
(4.15) 
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and 

E=^ + l; V = -^ (4.16) 
o o 

The Equation (4.14) is already known. It appears in the energy spectrum 
problem of an electron in a finite one-dimensional crystal in an uniform elec- 



tric field [j35|l , p6| . The electric potential Ex on the lattice gives a contribu- 
tion En on the lattice which appears in the left-hand side of Eq. (4.14). (We 
have divided by E the expression of the corresponding Hamiltonian) . For 
real D(q) the energy levels form the so-called Wannier-Stark ladder which 
has been investigated in great detail. Eq. (4.14) also appears in the problem 
of the diagonalization of the XX Hamiltonian in a Z field whose strength is 
position dependent |57| : 

TV 

H = ~ E + °M + i+ (h + RkK\ (4.17) 

k=-N 

where a x , a y and a z are Pauli matrices and h and R are constants. We will 



shortly present the solution of Eq. (4.14) as presented in Ref. |35|], it is given 
in terms of Lommel functions : 

<[ ILL 

^4 (-lf(m-n)! Tjfi + m-n) ( z ym+2n 

RmAz) = h n\{m-2n)\ T^ + n) U ) ^ 

which satisfy the identity 

RmA Z ) = (-l) m R m ,-^ m+ l(z) (4.19) 

The eigenvalues E of the equation (4.14) are given by the equation: 

Rl-i,i-e (2V(q)) = (4.20) 

There are (L — 1) solutions to the equation (4.20). To each of these solution 
corresponds an eigenfunction: 

F(n, q) = R n -i,i-E (2V(q)) F(l, q) ; (n = 2, . . . , L - 1) (4.21) 

Because of the identity (4.19) if E is an eigenvalue so is E where 

E = -E + L (4.22) 
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This implies that the eigenvalues E are distributed symmetric around |f. 
For L even, -| is an eigenvalue. For various properties of the zeroes of the 



Lommel polynomials see Ref. [SB]. To the eigenvalues obtained from Eqs. 
(4.16) and (4.20), one has to add the eigenvalue A = L5 obtained directly 
from Eq. (4.10) for q = 0. The general character of the spectrum for real 
V is the following. For not too large V, the eigenvalues in the center are 
equidistant, forming a simple ladder with unit spacing which results from 
the 1. h. s. of Eq. (4.14). At the upper and lower end of the spectrum 
the spacing increases. In this region the V-term in the equation becomes 
important. 

In the limit of large V or small 5 (i. e. for vanishing field in the electronic 
problem) one comes back to a simple tight-binding model with eigenvalues 

^ = 2cos^ (jfe = l,...L-l) (4.23) 

V Li 

and the corresponding eigenfunction is 

F(n,q) = —-^F(l,q) (4.24) 
sin 

This implies that for 5 = 0, the eigenvalues of the Hamiltonian are 

X(q) = 7 _2D(g)cos^ (Jfe = 1, 2, . . . L - 1; q = 0, 1, .L - 1) (4.25) 

to which we have to add A = 0, obtained directly from Eq. (4.10) in the 
case q = 0. From the expression (4.12) of D(q) we can obtain the conditions 
under which the spectrum is that of free fermions. The energy levels should 
correspond to two-fermionic excitations (see Sec. 3). If 



or _ Or 
a L a L 

one obtains 



(4.26) 



X(q) =l + G(q + k) + G(q-k) (4.27.a) 
where 

ll/2 ^ , : a R 



G(k) = - [(a R + a L )(a R + a L )] L/2 cos — + i^—^- sin — (4.27.b) 

V L a R + a L 
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which indeed looks like a two fermions excitation one of momentum (q + k), 
the other of momentum (q — k) with a total momentum q. The known 
examples (4.2), (4.4) and (4.6) satisfy the condition (4.26). Ha R /a L ^ a R /a L 
and if S = the spectrum is not given by free fermions, nevertheless the 
spectrum and wave-functions have a simple expression. 

We have checked on small chains that if the condition (4.26) is fulfilled the 
spectrum of the Hamiltonian (not only the one-hole subset) is indeed given 
by free fermions. We didn't look for the similarity transformation which 
would bring the Hamiltonian to the corresponding form. 

The problem of the open chain, also probably solvable, was not considered 
yet. We have only looked at small chains for the simple case where bn = b^ = 
c = 0. The levels are simply or doubly degenerate (the E = level is always 
doubly degenerate). This is in contrast to the free fermionic picture where 
all levels should be double degenerate (see Sec. 3). This is to be expected as 
long as the condition (4.26) is not satisfied. The surprise is that if one takes 



although the periodic chain is "fermionic" this is not the case for the open 
chain. 

Finally, we didn't touch the other sectors of the Hamiltonian (more than 
one hole). We will return to this problem in Sec. 7. 

5 The spectra in the continuum limit (closed 
ring) 

It is interesting to consider the continuum limit of the spectra. We first 
take the case 5 = 0. Then Eq. (4.14) becomes 




Or _ OR 

a L a L 



(4.28) 



( 



E 



2)F(n,q) 



F(n + 1; q) + F(n - 1; q) - 2F(n; q) 

d 2 F(n; q) 
dn 2 



V 



(5.1) 



with the boundary condition 



F(0;q)=F(L;q) = 



(5.2) 
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which implies 

E n 2 k 2 . , . 

--1 = -— (k=l,...) (5.3) 

and using Eq. (4.16): 



21,2 



A = 7 - 2D{q) + D(q)^- (k = 1, 2, . . . ; q = 0, . . .) (5.4) 



where (see Eq. (4.12)) 



D(q) = ^a R + a L )(a R + a L ) r , ^ ^ + ^ , ^ + ^ ; , „ 

(5.5) 

To Eq. (5.4) we have to add the value A = for q = (keep in mind 
Eq. (4.10)). The spectrum given by Eq. (5.4) is massive if 7 — 2D(0) > 0, 
massless if 7 = 2D(0) with a quadratic dispersion relation in the real part 
and a linear one in the imaginary part. 

We now consider the case 6^0. Using Eqs. (4.14) and (4.16), we get in 
the continuum: 

d 2 F(n,q) 



where 



with 



a + bn)F(n,q) (5-6) 

(5.7) 



dn 2 

7 — 5 — A 5 



D(q) ' D(q) 



F(0;q)=F(L;q) = (5.8) 
If Ai(z) and Bi(z) are two independent solutions of the Airy equation PU[ : 

— = zw (5.9) 

the general solution of Eq. (5.6) reads 

F(n, q) = C A Ai{a r 2/3 + b 1/3 n) + C B Bi{a r 2/3 + b 1/3 n) (5.10) 

where and Cb are arbitrary constants. 

The eigenvalues A (contained in the constant a) are determined from the 
boundary conditions (5.8). These equations have to be solved numerically. 
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In the infinite volume limit [L — > oo) the whole picture simplifies. Since the 
Airy function Bi(z) diverges for large z, we are left with the function Ai only. 
In order to satisfy the condition 

F(0,g) = (5.11) 

a b~ 2 / 3 has to be a zero of Ai(z). It is known |39[ that the zeroes of Ai(z) are 
all on the negative real axis in the complex z phase (they are tabulated in 



[j3"9fl), we give the first two values: C\ = —2.338 , c 2 = —4.087. This implies 
that 

A n = 7 - S - 2D(q)+ | c n | D\§ 6 2/3 (5.12) 

and thus the spectrum in the continuum is known. The eigenfunctions are 
given by Eq. (5.10) with Cb = and a fixed by Eq. (5.12). To the eigenvalues 
(5.12) one has to add for q = the eigenvalue A = that one gets from (4.10) 
for L infinity. At least in the domain where Or^, b^i, d^i and c are posi- 
tive the spectrum has to be massive. In this case one can use the stochastic 
interpretation of the Hamiltonian (3.10). The birth and pair-creation pro- 
cesses (see Eq. (2.4)) come with independent scales in the evolution operator 
determining a massive phase so that as long as 5 ^ one is in this phase. 
This statement is valid not only for the one-hole subset of the Hamiltonian 
but for all the subsets. 



6 The critical dynamics of the one-dimensional 
Ising model 

As an application of the calculations of spectra presented in the last section 
we consider the dynamics of the one-dimensional Ising model defined by the 
equilibrium probability distribution 

L 

P eq = Z- 1 ]J ( 6 .i) 

i=l 



where Z is the partition function: 



e n 



Si Si-\ 



Si=±l 



(6.2) 
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The two-point function has the following large distance behaviour 

< Si S i+R >~ e-M (6.3) 

where the correlation lenght ( and the mass \i are 

C -i = 2/ i , // = e"T (6.4) 

We can rewrite (6.2) in a sligthly different form: 

P e , = Z- l ne* ( - 1)( "'~ ,Hfl) =Z- 1 f[e^- 1 ^ (6.5) 

i=l i=l 

where = 0, 1. In the last equality of Eq. (6.5) we have performed a 

duality transformation. It was shown in Ref. f?j] that the following relations 
have to be satisfied by the rates appearing in the master equation (2.2) in 
order to have P eq as a steady state solution (detailed balance): 

Mo,o - r i,o = V> r o,i - Mil = r i,i - V> r o;o ( 6 - 6 ) 

-pO,0;a | T-\l,0;a . 2-nl,ha . T-\0,l;a rlA a la n\ 

where 

1 7,5 — 1 7,5 "+* 1 5,7 

p«,/3;a pCf:/3 p/8,c« 

7,5 7,5 5,7 

Let us assume that we give a set of rates compatible with Eqs. (6.6) - 
(6.7). This defines a Hamiltonian and the physical question one asks is the 
behaviour of the energy gap Eq in the thermodynamical limit as \i goes to 
zero (the Ising model in one dimension is critical at T = 0). The inverse of 
the energy gap gives the time correlation length. As shown in Ref. in 
several examples, two scenarios are possible. In the first, as \i — > 0, Eq stays 
finite. In this case we have no critical slowing down although we take local 
dynamics. In the second scenario, Eq vanishes like 

Eq ~ n z (6.8) 



with z = 2 (see also Refs. QT3 1 , |2~3" 1, |[24|| ) . It is our aim to check if for rates 



r"'f (3.10) given by the hole picture or spin-string rates f"'f given by Eq. 
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(3.25) one can obtain new solutions. Since in some cases (see Eq. (4.20) 
or (4.27)) the spectrum is complex it is interesting to see if the equilibrium 
state can be reached in an oscillatory mode. 

We first consider the holes picture. The Eqs. (6.6), (6.7) together with 
positivity give the following solution: 



b R = b L = c = ; a R = a L + rjfi 
d R = {J>{a L - if) ; d L = n(a L +rj) + i]fx 2 

where cll > and r] is a parameter satisfying the condition: 

a L > \ i] | 

Notice that 5 = (see Eq. (3.17)) and that 

a R _j_ d R 
a L d L 



(6.9) 
(6.10) 

(6.11) 



which implies that we do not have free fermions (see Eq. (4.26)) thus the 
model we consider is new. In order to compute the energy gap, we use Eqs. 
(5.4) and (5.5) to get: 

a L 2 



E, 



G 



7 



2j(a R + a L )(a R + a L ) = 



(6.12) 



where we have used Eqs. (3.17) and (6.9). As we see we get no oscillatory 
behaviour, the energy gap vanishes for [l — > again with z = 2 (see Eq. 
(6.8)) in agreement with the universality hypothesis. 

We now consider the rates f" ^ given by Eq. (3.25), in this case we find: 

8 ^ 2 A a u u c 

c = j- ^tOr , a R = a L , d R = d L = -b R = -b L = - 

l-/i 2 2 



(6.13) 



which gives: 



01 



r 



00 
10 



'--01 
00 



MO 
00 



P 10 
1 01 



- f 11 

- 1 01 
foi . 

1 10 > 



Ul 



S01 



- ^ 10 — 1 11 
plO i fiOl 



01 



+ r 



10 



_ f 10 _ Q 

- 1 00 — u 
foo i fill 
11 ' 00 



(6.14) 



Now we have free fermions and get back the Glauber model |fL2| . Since we 
want a R to be finite for /i — >■ 0, we choose: 

/i 2 A 



(A>0) 



(6.15) 
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where A is a constant and we get the known result 

E G = A/i 2 (6.16) 

Our results for Eq were obtained only within a sector of the Hamiltonian 
and one can always argue that other subsets of differential equations can 
give different results. This argument holds for Eq. (6.12) but not for Eq. 
(6.14) since in this case the Hamiltonian can be diagonalized completely . 



7 Conclusions 

We have shown how to bring a Hamiltonian describing a one-dimensional 
quantum chain with two states per site and depending on seven parameters 
to a block-diagonal form. The structure of the blocks is reminiscent of a 
Hamiltonian described by free fermions. The underlying algebraic structure 
of the problem is still a mystery for us. We have diagonalized the first 
block in the case of periodic boundary conditions. It turns out that Lommel 
polynomials play a fundamental role in the eigenvalues and eigenfunctions of 
the Hamiltonian. We have also clarified the continuum limit of the spectra. 
The case of the open chain was not considered but probably the method used 
in one special case can be extended to the general one. Hopefully the 



methods similar to the one used for the Glauber model [15] can be extended 
to all the blocks. If this is not possible, we will be left with an interesting 
example of partially integrable systems in condensed matter. 

We have paid special attention to various similarity transformations which 
connect stochastic observables not only in order to find properties of the spec- 
tra coming from positivity but also in order to understand how the method 
used in this paper can be extended to models with more than two states. 

In this paper we didn't do any specific calculations related to diffusion- 
reaction processes which would have implied solving also the nonhomogenous 
case of the differential equations since our purpose was to find the eigenvalues 
and eigenfunctions of the Hamiltonian. 

One basic question is if our Hamiltonian cannot be diagonalized using 
conventional Bethe-Ansatz methods. The answer is not obvious. We first 



have checked if one can't use the Baxterisation programme [ 4"0] . Before 
starting it we have checked if the matrices 

R k = H k + cl (7.17) 
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where H k is given by Eq. (3.10), c is an arbitrary constant and 1 the unit 
matrix, satisfy the braid group relations: 

RkRk+iRk = Rk+iRkRk+i (7-18) 

The answer is in general negative. If the answer would have been positive 
we would have had to look for some conditions to find an associative algebra 
for which the baxterization programme can be done. We have next tried 
the Reshetikhin criterium according to which the following relation has to 
be satisfied Q: 

[H k + H k+ x, [H k , H k+ i}} = X k — X k+ i (7.19) 

where X k is an arbitrary 4x4 matrix. The answer is again negative in general. 
One gets a positive answer only in what we called the free-fermionic case 
(the parameters satisfy the condition 5 = and Eq. (4.26)). As it is known 
the Reshetikin condition is very restrictive (the Hubbard Hamiltonian does 
not satisfy it) but the authors' knowledge about conventional integrability 
conditions stops here. 
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